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ABSTRACT 

Context. SiO maser emission occurs in the extended atmosphere of evolved stars and can be studied at high angular resolution. 
As compact, high brightness components they can be used as important tracers of the dynamics at distances close to the central 
star. The masers also serve as probes of the evolutionary path from spherically symmetric AGB stars to aspherical PNe. Very long 
baseline interferometry (VLBI) observations of Mira variables indicate that SiO masers are significantly linearly polarized with linear 
polarization fraction up to 100%. However, no information is available at high angular resolution for SiO masers in higher mass loss 
OH/IR stars. Theory indicates a different SiO pumping mechanism in higher mass loss evolved stars. 

Aims. We extend the VLBI SiO maser studies to OH/IR stars. The observations enable us to understand the SiO pumping mechanisms 
in higher mass loss evolved objects and compare those with Mira variables. Additionally, polarimetric observations of SiO masers 
help us to understand the magnetic field strength and morphology and to distinguish between conflicting polarization theories. 
Methods. The 43 GHz SiO maser observations of the OH/IR star OH 44.8-2.3 were performed with the VLBA in full polarization 
spectral line mode. Auxiliary EVLA observations were performed to allow for the absolute calibration of the polarization angle. 
The Zeeman splitting was measured by cross correlating the right and left circular polarization spectra as well as the S-curve fitting. 
Additionally, we analyzed the 1612 MHz OH maser observations of OH 44.8-2.3 from the VLA archive. 

Results. The SiO masers of OH 44.8-2.2 form a ring located at ~5.4 AU around the star. The masers appear to be highly linearly 
polarized with fractional linear polarization up to 100%. The linear polarization vectors are consistent with a dipole field morphology 
in this star. We report a tentative detection of circular polarization of ~0.7% for the brightest maser feature. The magnetic field 
measured for this feature corresponds to 1.5+0.3 G. Additionally, the distribution of the 1612 MHz OH maser emission could indicate 
an elongated morphology. 

Key words. Stars — magnetic fields — polarization — masers 



1. Introduction 

SiO maser emission associated with Asymptotic Giant Branch 
(AGB) stars occurs at the inner region of Circumstellar 
Envelopes (CSEs). Maser emission has been detected in several 
vibr ationally-excited rotational transitions of the SiO molecule 
(e.gJHumphrevs et aljfl997t iPardo et"al]|1998tlSoria-Ruiz et al] 
2001. A series of very long baseline array (VLBA) high resolu- 
tion observations of the SiO masers around Mira variables have 
shown that the maser emission is confined to a region, sometimes 
ring-shaped, be tween the stellar photosphere and the dust forma - 
tion zone (e.g. ICotton et alJl2008l 120061: iDiamond et ail 1 19941) . 
On scales probed by very long baseline interferometry (VLBI), 
the masers are confined to localized spots with lifetime of a few 
months. The resolution of the VLBI observations of SiO masers 
is a small fraction of the ring. The above points imply that SiO 
masers can be used as powerful probes of the processes that drive 
the mass loss and dynamics in the inner region of the CSEs. 

SiO masers can also be important for understanding the evo- 
lutionary path from the spherically symmetric AGB stars to as- 
pherical Planetary Nebulae (PNe). Magnetic fields likely play 
an important role in shaping the CSEs of evolved stars (e.g. 
iGarcia-Segural [T997h . Polarization observations of circumstel- 



lar masers enable us to determine the magnetic field strength 
and morphology at different distances from the central stars. 
Observations indicate that SiO masers are significantly linearly 
polarized, with polarization vectors mainly tangential to the 
maser ring. For example, multi-epoch polarization VLBA obser- 
vations of the Mira star TX Cam revealed the linear polarization 
morpholo gy tangential to the p rojected shell of SiO maser emis- 
sion (e.g. IKemball et al J 120091) . However, Cotton et al. (2008), 
observed five Mira variables with the VLBA at 43 GHz, but not 
all sources showed a tangential linear polarization morphology. 
It remains therefore unclear whether the tangential linear polar- 
ization morphology is a generic property of the SiO maser region 
of evolved stars. Th e circular polarization of the masers is in th e 
range 3%-5% (e.g. iBarvainis et alJ[l987l: IKemball eTall 12009). 
Polarimetric observations of a sample of ev olved stars indicate 
an average magnetic field of several gauss dHerpin et al.| [2006: 
Kemball & Dia mondlfT997h . The SiO maser magnetic fields are 
consistent with those measured using H2O and OH masers fur- 
ther out in the envelope as suming a dipole (B oc r~ 3 ) or solar-type 
(B oc r~ 2 ) magnetic field dVlemmings et alJl2005l) . 

VLBI studies of SiO masers have so far focused on Mira 
variables and a handful of post- AGB objects and water fountain 
sources (e.g. lCotton et al.ll2008t IPesmurs et al1l2007tllmai et all 
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I2005h . and as far as we can tell no high resolution VLBI ob- 
servation has been performed to image SiO masers in OH/IR 
stars. These objects have larger CSEs and much longer pe- 
riods up to 2000 days co mpared to those of Mira variables 
dHerman&Habindll985bh The v are strong 1612 MHz OH 
maser emitters dBaud et al] [T979). The stars are surrounded by 
thick dust shells which makes them optically obscured. The SiO 
maser pump mechanism in OH/I R stars may be d ifferent from 
that operating in Mira variables dDoel et alj|1995l) . Due to the 
higher mass-loss of these objects the CSE is denser at least an 
order of magnitude more than those of Mira variables. 

Here, we report the SiO maser polarimetric observations of 
the OH/IR star OH44.8-2.3 with the VLBA. The observations 
enable us to obtain the spatial distribution of the SiO maser fea- 
tures in OH/IR stars for the first time. Additionally, our experi- 
ment probes the magnetic field strength and morphology in the 
SiO maser region of OH/IR stars and compares them with those 
of Mira variables. We also reduced the 1612 MHz OH maser 
observations of this star from the VLA archive in order to deter- 
mine whether there exists any large scale asymmetry in the OH 
maser shell of this star. 

OH 44.8-2.3 is identified as a moderate OH/ IR star with a 
warm and thin CSE and a period of 534 days dGroenewegenl 
11994 . The mass loss of this star is estimated to be ~4.6xl0 
M Q /yr dde Beck et alj2010h . The OH masers of this source show 
a double peak profile at -88.9 and -53.8 km s _I with a peak 
flux density of 14.4 Jy for the re d shifted part of the spectrum 
dEngels & Jimenez^ Esteban 2007). A distance of 1.13±0.34kpc 
was measured for this star using the phase-lag method based 
on the characteristic 16 12 MHz double peak profile of the OH 
masers of this source dVan Langevelde et al.l [19 90). The H2O 
masers of this star were detected by En gels et alJ J 19861) at ~5 Jy. 
However, further single dish monitoring of this source did no t 
result in any detection dShintani et alj|2008t Ikim et alj |2010). 
Th e v=l, J — 1 — > SiO masers of this star wer e observed 
bvlNvman et all d 19981) with a flux density of 17.6 Jy. Ki m et alJ 
(2010) observed the v= 1,2,7 = 1 — > emission from this star 
with an integrated flux density of ~25 Jy and ~28 Jy, respec- 
tively. 

The outline of the paper is as follows: the observations are 
described in Sec. [2] Before giving the results in Sec. [4] we intro- 
duce the necessary background on maser polarization in Sec. [3] 
We discuss the interpretation of the results in Sec. [5] This is fol- 
lowed by conclusion in Sec . [6] 



2. Observations 

We observed the v=l, J = 1 — > SiO maser emission to- 
ward OH 44.8-2.3 on 6 July 2010 using the NRAO Q Very 
Long Baseline Array (VLBA) operating in the 43 GHz band. 
Additionally auxiliary Extended Very Large Array (EVLA) in- 
terferometer observations were performed on 2 July 2010 and 
1 1 July 2010 to allow for the absolute calibration of the electric 
vector polarization angle (EVPA). Furthermore, we retrieved the 
1612 MHz OH maser very large array (VLA) observations of 
OH 44.8-2.3 from the NRAO archive. The VLBA, EVLA and 
VLA observations are discussed below: 



1 The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



2.1. VLBA observations and reduction 

The data were recorded in dual circular polarization spectral line 
mode, which generates all four polarization combinations in the 
correlator. The DiFX correlator was used with a bandwidth of 4 
MHz and 1024 spectral channels, which results in 0.03 km s _1 
spectral resolution. The observations were performed in one 
spectral window centered at a fixed frequency corresponding to 
the v=l, J — 1 — » SiO maser transition at a rest frequency of 
43. 12208 GHz and a stellar velocity of -72 km s _1 with respect to 
the local standard of rest (LSR), determined from OH maser ob- 
servations. The total observing time was 6 hr, balanced between 
the target source, OH 44.8-2.3, and the continuum calibrators, 
J2253+1608, J1751+0939 and J 1800+3848. We achieved a spa- 
tial resolution of 0.5x0.2 mas. 

We used the Astronomical Image Processing Software 
Package (AIPS) to perform the rest of the calibration, editing and 
imaging of the data. The first calibration steps were performed 
on the data set with modest spectral resolution (128 channels). 
The solutions were then applied to the high spectral resolution 
data set (1024 channels). Parallactic angle corrections were per- 
formed on all calibrators. The North Liberty (NL) and Saint 
Croix (SC) antennas had to be flagged due to the bad weather 
conditions at the time of the observations. 

Since the SiO molecule is non-paramagnetic, the fractional 
circular polarization will be small (m c ~ 1% — 3%). In order 
to preserve the low stokes V signature, we performed ampli- 
tude calibration of OH 44.8-2.3 using the AIPS task 'ACFIT'. 
In this method, each circular polarization autocorrelation spec- 
tru m was calibrated independently using the method described 
bv lReid et al.l d!980t) . A template auto correlation spectrum was 
selected for ~20 min interval from the Los Alamos antenna at a 
time range with sufficiently high elevation. The template spec- 
trum was fitted to all other total power spectra from all antennas 
and the relative gains of the antennas as a function of time were 
determined. We used the system temperature measurements pro- 
vided with the data to perform the amplitude calibration for the 
calibrators. 

The complex band pass solutions were obtained using 
J2253+1608 on each data set separately. Fringe fitting for the 
residual relay and rate for the parallel hand data was done using 
the continuum calibrators J2253+1608 and J1751+0939. The 
fringe rate solutions for OH 44.8-2.3 were determined on the 
strongest feature. The delay offset between the right and left cir- 
cular polarizations at the reference antenna was determined for 
J2253+1608 on the Los Alamos - Fort Davis baseline. The so- 
lutions were transferred to all other baselines and sources. The 
polarization leakage terms were determined using J2253+1608. 
After the solutions were applied to both data sets image cubes 
could be created with 50 yuarcsec pixel spacing. The rms noise 
in the modest spectral resolution and the high spectral resolution 
individual channel images corresponds to 8 mJy and 26 mJy, re- 
spectively. 



2.2. EVLA observations and reductions 

The absolute phase difference between the right and left circu- 
lar polarizations at the reference antenna is not known, which 
implies that the absolute EVPA of linearly polarized emis- 
sion cannot be measured accurately in our VLBA observa- 
tions. Therefore, auxiliary EVLA observations were performed 
to measure the absolute polarization angle of the polarization 
calibrators with respect to a primary polarization calibrator with 
known polarization angle. 
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The continuum sources J2253+1608 and J1751+0939 were 
observed as transfer calibrators on 2 July 2010 and 11 July 
2010, respectively. The observations were performed in contin- 
uum mode and full polarization, using two 128 MHz spectral 
windows. The data reduction was performed for each spectral 
window separately. We achieved a beam size of 2.6x1.7 arcsec 
during the observations. 

We reduced the data according the standard EVLA data cali- 
bration and imaging recipe in the appendix of AIPS Cook Book. 
Since the delays are not set accurately in EVLA observations, we 
solved for the phase slope using the AIPS task FRING. A ~ 1 min 
solution interval was chosen on J2253+1608 and J1751+0939 
to solve for the phase slope as a function of frequency. The 
bandpass solutions were determined using J2253+1608 and 
J1751+0939. The flux density of the primary flux calibrators 
3C286 and 3C48 was determined using the model provided 
within the AIPS software. The AIPS task 'CALIB' was used 
to perform phase calibration for 3C286 and 3C48 using a model 
of the sources. Subsequently, the phase solution was determined 
for the secondary calibrators J2253+1608 and J1751+0939. The 
amplitude of the secondary calibrators was adjusted with respect 
to the flux density of the primary flux calibrator. We used 3C286 
and 3C48 for polarization calibration. The delay difference be- 
tween the right and left circular polarizations were determined 
using the AIPS task 'RLDLY'. The feed parameters for each po- 
larization were determined both in continuum and spectral mode 
since the EVLA feeds have significant variations in frequency. 
Images were produced in each spectral window in Stokes Q, 
U, I adopting a pixel spacing of 0.2 arcsec. The rms noise for 
J2253+16 and J175 1+9839 were ~8 mJy and ~1 mJy, respec- 
tively. 

The absolute EVPA of J225 1608 and J 175 1 +0939 were mea- 
sured with respect to 3C48 and 3C286, respectively. The fi- 
nal absolute EVPA was calculated as the mean of the EVPA 
in each spectral window. This corresponds to the measured 
EVPA of 62° and -44° for J2253+1608 and J 175 1+0939, respec- 
tively. Fig.Q]shows the EVPA of the linearly polarized emission 
for J2253+1608 for both EVLA and VLBA observations. The 
EVPA in the VLBA image is rotated by 15°, which corresponds 
to the difference in the EVPA between the EVLA and VLBA 
images. We note that the linearly polarized emission of the other 
polarization calibrator, J 175 1+0939 shows two polarized com- 
ponents separated by ~90° in the VLBA image. Therefore we 
were not able to use this source as a robust EVPA calibrator. 



2.3. VLA observations of the 1612 MHz OH masers of 
OH 44.8-2.3 

We found the previous observations of the 1612 MHz OH masers 
of OH 44.8-2.3 observed on 25 September 1983 using the Very 
Large Array (VLA) from the NRAO archive. The observations 
were performed under the project name 'AH 127' in the A con- 
figuration with the largest spacing of 36.5 km, which gives a 
resolution of 1". The band width of 1.5 MHz was used for the 
observations which gives a velocity resolution of 1.4 km s _1 at 
1612 MHz. 3C286 was observed as a flux density calibrator and 
1741-038 was used as the secondary phase calibrator. The data 
reduction was performed following the recipe provided in the 
AIPS cookbook. 

Fig. displays the OH maser map of OH 44.8-2.3. The 
emission covers a velocity range of -87 to -56 km s _1 . The fig- 
ure shows the blue- and red-shifted peaks as well as emission 
summed over several channels close to the stellar velocity. 



3. SiO maser polarization theory 

The SiO maser emission at 43.12208 GHz involves the rota- 
tional transition / = 1 — > in the excited vibrational state v=l. 
Since the SiO molecule is non-paramagnetic, its response to the 
magnetic field is weak. This implies that the Zeeman splitting is 
smaller than the line width. 

3. 1 . Linear polarization 

SiO masers present high f ractional linear polarization 
(e.g. iKemball & Dia mond 199J). It was shown that anisotropic 
background radiation from the central star generates anisotropic 
pumping, potentially producing a high de gree of linear po- 
larization in the SiO emission region (e.g. IWatsonl l2002l) . In 
the case where the SiO masers are radiatively pumped, the 
magnetic sub states may be anisotropically populated, which 
can ge nerate high linear polariz ation fraction for saturated 
masers [We stern & Watson (1983). However, the results from 
iNedoluha & Watsorj (1 19901) show that the polarization vectors 
still trace the direction of the magnetic field despite the fact that 
the linear polarization may originate largely from anisotropic 
pumping. 

3.2. Potential non-Zeeman effects for circular polarization 

The non-Zeeman mechanism prohibits one to interpret circular 
polarization as a measure of the magnetic field strength. This 
stems from the competition between the stimulated emission rate 
(R), the Zeeman coefficient rate (gQ) and the radiative decay rate 
(r). The circular polarization can be created by the change of 
the axis of symmetry for the molec ular quantum states when th e 
condition gQ R > F is satisfied (INedoluha & Watsodfl994T) . 
For some distance along the maser path gQ > R is satisfied 
and the magnetic field is the quantization axis. As the radiation 
propagates farther into the maser path and the rate of stimulated 
emission becomes larger, the inequality R > gQ is satisfied and 
the molecule interacts more strongly with the radiation. This 
implies that the axis of the symmetry of the molecule changes 
from parallel to the magnetic field to parallel to the direction of 
propagation. The resulting circular polarization would have the 
anti-symmetric profile which resembles to that produced by the 
ordinary Zeeman effect. The intensity dependent circular polar- 
ization can be higher than that created by the Zeeman effect by 
a factor as large as 1000. This scenario is an inherent part of 
the radiative transfer process and does not require anisotropic 
pumping or high fractional linear polarization. 

[Wiebe & Watson) (1 19981) introduce yet another non-Zeeman 
effect in which the propagation of a strong linear polarization can 
create circular polarization if the condition gQ » R > T is satis- 
fied. The circular polarization can be generated if the magnetic 
field orientation changes along the maser propagation direction. 
The circular polarization produced from this scenario is on av- 
erage ~ -j- where nil indicates the linear polarization fraction. 
For a typical linear polarization fraction around 30%, the circu- 
lar polarization results from the formula above corresponds to 
~2-3% whi ch is in agreement wit h the observed circular polar- 
ization (e.g. lBarvainis et aT] [T987). however in individual maser 
features the circular polarization fraction can go up to ~20%. 
This mechanism would then generate circular polarization with 
only 10 - 20 mG fields. Moreover, if there is no velocity gradi- 
ent along the maser path, the circular polarization profile would 
have the shape of the S-curve. However, the intrinsic profiles can 
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Fig. 1. Linearly polarized map of the polarization calibrator J2253+16 obtained with the EVLA (left panel) and the VLBA (right 
panel). The vectors are rotated by 15° in the VLBA image. 




Fig. 2. The 1612 MHz OH maser map of OH 44.8-2.3 obtained with the VLA. The top-left and bottom-right panels show the blue- 
and red-shifted peaks. The red-shifted peak exhibits the highest flux density of 8. 1 Jy. The contour levels are at 1 .62, 3.24, 4.86, 6.48 
and 8.1 Jy. The top-right and bottom-left panels exhibit emission summed over several channels close to the stellar velocity. The 
contour levels are at 0.005, 0.01 1, 0.016, 0.022 and 0.027 Jy. The star symbol indicates the position of the peak in the red-shifted 
emission which likely indicates the position of the central star. 
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be distorted as a result of saturation, together with velocity and 
magnetic gradients along the amplification path. A correlation 
between the linear and circular polarization fraction is expected 
for the individual features. 



4. Results 

4. 1 . Total intensity 

Fig.[3]shows the total intensity spectrum of the SiO maser emis- 
sion towards OH 44.8-2.3. The emission covers a velocity range 
of ~1 km s _1 . The SiO emission actually implies a stellar ve- 
locity 3-4 km s offset from the stellar velocity of -72 km s 
based on OH maser observations of this star. Fig.|4]displays the 
SiO maser emission map of OH 44.8-2.3 summed over all ve- 
locity channels covering emission. The maser features form a 
partial ring of 4.75 mas corresponding to 5.4 AU assuming a 
distance of 1.13 kpc dVan Langevelde et al.ll 1990T> . The masers 
appear to be absent from the North-East and South- West part of 
the ring. The ring pattern observed in the SiO maser region of 
this star implies that the masers are tangentially amplified. This 
is typically interpreted to indicate that large velocity gradients 
exist in the SiO maser region which prohibits radial amplifica- 
tion. The peak flux density for the maser features are displayed 
in TableQ] Feature 1 exhibits the largest flux density of 2.8 Jy. 



4.2. Linear polarization 

The polarization morphology is shown in Fig. In this plot, 
the polarized emission is plotted as vectors with a length pro- 
portional to the polarization intensity. The position angle of the 
vectors corresponds to the EVPA of the emission. All stokes pa- 
rameters are summed over frequency before making this plot. 
The maps of stokes I, Q, and U are available via CDS. The 
background contours represent the total intensity image. The 
fractional linear polarization for the individual maser features 
is shown in Table [TJ The average linear polarization fraction for 
the SiO maser features of this star is ~30%. In particular for fea- 
ture 8 it reaches 100%. This likely implies that the (radiative or 
collisional) pumping of SiO masers is not isotropic. 

4.3. Circular polarization 

The instrumental gain differences between the RCP and LCP 
spectra, can cause a residual scaled down version of the Stokes 
I spectrum to appear in the circular polarization profile. The 
Zeeman splitting for the stokes V spectra of the SiO maser fea- 
tures of OH 44.8-2.3 can be determined by measuring the resid- 
ual curve having the shape of the d erivative of the total inte nsity 
spectrum in the least square sense (Trol and & Heil es 1982): 



dl 

V(v) = a x — + b x /(v) 
dv 



(1) 



Where a is the measure of the circular polarization due to the 
Zeeman splitting and b denotes the instrumental differential am- 
plitude. This method assumes that intrinsically equal flux is 
present in the RCP and LCP spectra, which implies that the cir- 
cular polarization spectra have the antisymmetric S-curve pat- 
tern. However, theoretical models of circular polarization pro- 
duce asymmetric profiles de pending on the velocity gradient 
along the amplification path (Wiebe & Watson] fl998). We find 
a tentative detection of the circular polarization for the brightest 
maser feature in Fig. [4] (feature 1 from Table[TJ. Fig.|6]displays 



the circular polarization together with the total intensity spectra 
for this feature. The figure also shows the fit to the circular polar- 
ization profile after removing the scaled down replica of stokes 
I. The circular polarization fraction is measured as: 



V m 



x 100 



(2) 



Where V„ MV corresponds to the maximum of the fit to the ob- 
served stokes V spectrum (Eq.Q] a x ^). For the modest spec- 
tral resolution data (128 channels), we measure circular polar- 
ization of m c ~ 0.7 + 0.2%. However, due to the increased 
noise in individual channels in the high spectral resolution data, 
we can not confirm the detection. The magnetic field derived 
from circular polarization co rresponds to the following equation 
(Kemball & Diamon d! 19971 derived from Elitzur (1996)): 



B = 3.2 x m c x Av D x cos 6 



(3) 



Where m c , Avq and 9 indicate the fractional circular polariza- 
tion, the maser line width and the angle between the magnetic 
field and line of sight, respectively. The full width half maxi- 
mum line with for the stokes I spectrum of feature 1 corresponds 
to ~ 0.8 km s _1 . Using the preceding relation a magnetic field of 
1.8+0.5 G is derived for feature 1. 



4.3.1. Cross-correlation method and magnetic field strength 

Alt ernatively, we use th e cross-correlation method introduced 
by [Modiaz et al. (2005) to measure the magnetic field due to 
the Zeeman splitting. In this method the right circular polariza- 
tion (RCP) and the left circular polarization (LCP) spectra are 
cross-correlated to determine the velocity splitting. The mag- 
netic field is determined by applying the Zeeman splitting co- 
efficient for SiO masers. This method can even work for com- 
plex spectra, assuming that the velocity shift is the same over 
the spectrum; which means the magnetic field strength and di- 
rection is constant in the maser region. The sensitivity of this 
method is comparable to the S-curve method, where the stokes V 
spectrum is directly used for measuring the mag netic field (e.g. 



IVlemmings et alfcOOUlFiebig & GuesterJl9f 



iagr 



Th e Zeeman rate for SiO is given by iNedoluha & Watson! 
(1990) as: 



§-*»(§) 



(4) 



The equation above implies the Zeeman splitting coefficient 
for the J = 1 — > 0, v=l transition of the SiO maser as 
8.7xl0~ 3 km s _! G We present the result of the cross- 
correlation method in Table Q] We measured a magnetic field 
of 1 .5+0.3 G for feature 1 . However, for the other features of the 
SiO emission of this star we only place upper limits at 3<x level. 

4.4. OH maser observations of OH 44.8-2.3 

The VLA results on OH 1612 MHz emission from OH 44.8- 
2.3 were not presented directly in the scientific lit erature before, 
althou gh the original distance d etermination in Herman et al.l 
(1986) and the one we use from IVan Langevelde et alJ d 19901) 
depend directly on it. We have access to the results of the origi- 
nal analysis in which a shell radius of 1.302" + 0.023" was ob- 
tained (van Langevelde, private communication). The new pro- 
cessing we performed gave consistent images with these previ- 
ous results, but we note that the OH at intermediate velocities 
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Velocity (km/s) 

Fig.3. The total intensity spectrum of OH 44.8-2.3 for the v=l, J — 1 — > SiO maser emission obtained with the VLBA. 
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RA 


Dec 


V 


I 


m, 


B 




19 21 


09 27 


km s -1 


Jy 


% 


G 


1 


36.63700337 


56.5199500 


-69.4 


2.9 


8 


1.5 ±0.3 


2 


36.63710476 


56.5203000 


-69.8 


0.7 


26 


<3 


3 


36.63722303 


56.5180500 


-67.8 


0.5 


23 


<9 


4 


36.63725683 


56.5169000 


-65.9 


0.9 


36 


<3 


5 


36.63728386 


56.5161500 


-67.7 


1.2 


34 


< 2 


6 


36.63731427 


56.5152000 


-68.3 


0.7 


48 


<2.5 


7 


36.63727710 


56.5141000 


-66.6 


0.3 


50 


<3 


8 


36.63717234 


56.5091500 


-70 


0.6 


100 


<2.7 


9 


36.63715545 


56.5113500 


-68.7 


0.4 


50 


< 12 


10 


36.63706420 


56.5106000 


-71.2 


0.6 


42 


<3 


11 


36.63698310 


56.5111500 


-67.2 


0.6 


18 


<36 


12 


36.63680399 


56.5109500 


-67.7 


0.5 


30 


<6 


13 


36.63673641 


56.5116000 


-67.9 


0.3 


69 


< 7 


14 


36.63661475 


56.5156000 


-69.0 


0.7 


46 


<2.5 


15 


36.63682427 


56.5214500 


-69.2 


0.24 




< 12 



Table 1. Results of the Magnetic fields determination of the SiO maser features of OH 44.8-2.3 obtained from the VLBA observa- 
tions. 



is quite elongated (Fig. 0. The inferred size of the region is of 
course much larger than the SiO, spanning ~1471 AU compared 
to 5.4 AU. One should note that with typical outflow velocities 
of ~16 km s _1 and 5 km s _1 for OH and SiO masers, the dynam- 
ical timescales are also much different: ~440 years for the OH 
and ~5 years for the SiO. 

5. Discussion 

5. 1 . Linear polarization 

In order to the compare the two conflicting polarization theo- 
ries (a magnetic versus anisotropic pumping origin of the linear 
polarization implying a Zeeman or non-Zeeman origin of the 
circular polarization), we need to determine in which theoretical 
regime the masers were observed. This includes the determina- 



tion of the stimulated emission rate (R), the Zeeman rate (gQ), 
and the collisional and radiative decay rate F. The radiative de- 
cay rate for the v= 1 , J — 1 — » transition of the SiO molecul e 
is estimated to be 5 s" 1 dKwan & ScoviUd[l974T: [Elitzurll 19921) 
The s timulated emission rate (R) is estimated by (Kemball et al. 
2009). 



R = 23 ( TB m ) 
\2x 10 10 K/ 



da 

10- 2 sr 



(5) 



Where Tg is the maser brightness temperature and dQ. is the 
estimated maser beaming angle. We obtained an estimate for 
the brightness temperature for feature 1 (Fig. [H> as 5xl0 10 K. 
We ad apt an estimated beaming angle of 10~ 2 sr (Kembal i et al.l 
2009). Replacing these values in eq. [5] results in R=57.5 s -1 . 
Substituting the value of ~1.5 G obtained for feature 1 (TableQ]) 
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Fig. 4. The VLBA map of v=l, J — I SiO maser emission towards the OH/IR star OH 44.8-2.3. The features are color coded 
according to the flux density (Jy / beam) integrated over all velocity channels. 



in eq. H implies g£l -1900 s" 1 . Therefore, we are in a regime 
where gQ » R > F. In this regime the linear polarization vec- 
tors appear either parallel or perpendicular to the projected mag- 
netic field, depending on the angle between the ma gnetic field 
direction and the line of sight (lGoldreichetalJl973l) . Taking the 
EVPA of -50° for feature 8 which has the highest linear polar- 
ization fraction, this implies that the magnetic field direction is 
either parallel (-50°) or perpendicular (40°) to the linear polar- 
ization vectors. In either cases, the EVPA vectors likely indicate 
a large scale magnetic field in the SiO maser region of this star. 



The linear polarization morphology of this star is not con- 
sistent with the dominant tangential polarization morpholo- 
gies that were repor ted for other evolved stars (e.g. TX Cam; 
Kem ball et al. 2009). Based on p revious observations of the 
SiO maser region o f TX Cam dKemball & Diamond! 119971; 
iDesmurs et alJl2000l) . they conclude that the tangential polar- 
ization morphology could be an inter cycle property of the SiO 
maser emission toward this star. Instead of tangential polariza- 
tion morphology, the EVPA vectors could indicate a bipolar 
magnetic field morphology for the SiO maser region of OH 44.8- 
2.3 (Fig. [5]). However, we can not rule out the other, more com- 
plex, morphologies. Polarimetric observations of the OH and 
H2O masers of this star are required to clarify the magnetic field 
morphology of this star. A dipole field morphology is reported 
for the supergiant VX Sgr usi ng the polarimetric obs ervations of 
high frequency SiO masers dVlemmings et al.ll201 lb . This mor- 
phology was consistent with the the dipole magnetic field in- 
ferred from H2O and OH maser regions at much larger distances 
from the central star ( Szvmczak & Cohenll99"7t lSzvmcza k et al] 



l200lHVlemmings et alJ^OOir 



5.1 .1 . The effects of anisotropic pumping 

The SiO maser features of this star exhibit high fractional 
linear polarization up to 100% (Table [TJ, which likely indi- 
cates that the pumping of the masers is highly anisotropic. 
Such high linear polarization fractions were rep orted previously 
(iKemball & Diamond! 1 9971: ICotton et al.ll2006f) . The masers oc- 
cur close to the central star and the infrared radiation incident 
on the masering region is highly anisotropic. This implies that 
the magnetic sub states are pumped unequally. Additionally, the 
medium in which they are located likely has large velocity gra- 
dients, wh ich leads to p r eferred direction fo r the p ropagation 
of masers (IWatsonll2002l) . I Western & Watsonl (|l983) show that 
even for collisionally pumped SiO masers, the magnetic sub 
states are anisotropically pumped and high linear polarization 
would be expected for saturated and unsaturated masers. In the 
absence of magnetic fields, collisional pumping causes linear po- 
larization along the radial direction. In contrast, radiative pump- 
ing would produce linear polarization vectors tangential to the 
SiO maser ring. However, our observations show that magnetic 
fields are present in the SiO maser ring of OH 44.8-2.3, which 
implies that both magnetic field and anisotropic pumping influ- 
ence the linear polarization morphology. Therefore, we can not 
distinguish between the radiative or collisional pumping from 
the observations. 



5.1.2. Maser saturation 

The degree of saturation is the ratio of the rate R for stimulated 
emission to the loss rate F. From the observations we measured 
the stimulated emission rate for feature one as R~57.5 s _1 which 
implies ^ ~ 10. This indicates that the medium is marginally 
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Fig. 5. Contour plot of the Stokes I image at levels [1, 2, 5, 10, 20, 40, 80, 100]% of the peak. Vectors are overlaid proportional to the 
linearly polarized intensity (on a scale 1 mas = 1.25 Jy beam -1 ) and drawn at a position angle of the EVPA. All Stokes parameters 
(I,Q,U) are summed over velocity. 
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Fig. 6. Total power (I) and circular polarization (V ) spectra of the brightest SiO maser feature of OH 44.8-2.3. The dashed line is 
the fit to the observed V-spectrum . The V spectrum is shown after removing the scaled down replica of stokes I. 

saturated. In contrast, a saturation of f ~ 30 is required to ( 1994) show that the low level circular polarization observed for 
achie ve 70% linear polarizati on fraction in the non-Zeeman see- SiO masers can be due to the saturation of masers. 
nario dWatson & Wvldll2001[). In the stand ard Zeeman interpre- 
tation scenario proposed by lElitzun (119961) the polarization solu- 
tion does not depend on the saturation level . iNedoluha & Watsonl 
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5.1.3. Jet-like features 

Even though the majority of the SiO maser features of OH 44.8- 
2.3 are coherent and fit the ring model, features 8 and 15 do not 
conform with the ring pattern. A striking characteristic of fea- 
ture 8 is the 100% linear polarization fraction. Additionally, the 
linear polarization vectors of this feature are in the perpendicu- 
lar direction of the overall l inear p olarization morphology of this 
star (Fig. l5l). ICotton et al.l d2006l) report structures aligned with 
the direction of the photosphere with magnetic field morphology 
along the features. They refer to them as jet-like features. They 
explain that these features likely form due to the masering region 
being dragged along the magnetic field. For the case of OH 44.8- 
2.3, the 100% linear polarization fraction could imply that this 
feature is a dynamic part of the envelope. Observations of this 
feature close in time to study the proper motion is necessary to 
understand the nature of this feature. 

5.2. Circular polarization 

We report a tentative detection of ~0.7% for the circular po- 
larization fraction for feature 1 (Table The measured mag- 
netic field corresponds to 1.5+0.3 G for this feature. We note 
that the magnetic field measured from the circular polarization 
fraction (Sec. I4.3U .8±0.5 G) and the one measured from the 
cross-correlation method (Sec. I4.3.T1 1.5+0.3 G) are consistent 
within errors. Therefore the correspondence of the field mea- 
sured from the cross-correlation method with the field determi- 
nation through the circular polarization fraction confirms that the 
cross-correlation method works properly.The non-Zeeman effect 
due to intensity dep e ndent circular polarization introduced by 
iNedoluha & Watson! d 19941) is ruled out since gQ » R in our 
observations. However, since we are in a regi me where gCl » 
R > r the non-Zeeman effect introduced by IWiebe & Watsonl 
( 1998) is applicable. For the 7% linear polarization fraction mea- 
sured for feature 1 (Table [TJ, the generated circular polarization 
due to this effect is 0. 12%. This implies that the measured circu- 
lar polarization for feature one is about 6 times high er than the 
estima ted value from the non-Zeeman effect. IWiebe & Watsonl 
(1998) show that if the circular polarization is higher than the 

ffl: 

average of -4-, the circular polarization stems from other effects, 
probably due to the Zeeman effect. Therefore, it is likely that 
the circular polariz ation of this star originat es from the Zeeman 
splitting. However. IWiebe & Watson (1998) explain that the av- 
erage circular polarization in individual features can go up to 
20%. Additionally, the determination of Zeeman or non-Zeeman 
effect from the observed circular polarization profile remains in- 
conclusive, since both models predict the similar anti-symmetric 
S -curve profile. 

5.3. CSE morphology and magnetic field 

The SiO maser features of OH 44.8-2.3 exhibit two opposite 
arcs (Fig. H]i. Additionally, the OH masers of this star show that 
the masers have an elongated shell morphology in the direction 
where there is a gap in the SiO maser emission (Fig. [2]i . We 
note that the OH masers occur on much larger scale (~1471 AU) 
than the SiO masers (5.4 AU) around the star. It is therefore not 
obvious that both deviations from symmetry are related, but if 
they are, this is an indication that there is a mechanism at work 
that can support the asymmetry on many scales. Similar gaps or 
opposite arcs were seen recently in monito ring observatio ns of 
the SiO masers of R Cas with the VLBA dAssaf et alJl201lh . Of 



course, only with one epoch of observations we can not deter- 
mine whether the SiO maser arcs in OH 44.8-2.3 are consistent 
in time since the region where SiO masers occur is expected to 
be highly dynamic. Therefore, multi-epoch observations of the 
SiO masers of this star is essential to probe the motion of the 
SiO maser features. 

As mentioned earlier, the magnetic field morphology for 
the SiO maser region of this star is either parallel or perpen- 
dicular to the linear polarization vectors. Interestingly, the di- 
rection of the magnetic field is parallel or perpendicular to 
the location of the gaps in the SiO maser ring and the OH 
maser extent. There thus appears to be a global preferred direc- 
tion of the outflow imposed by the magnetic field in the CSE 
of this star. Such asymmetric outflows have already been re- 
ported for the CSE of evolved stars. For example, collimated 
H2O maser jets have been observed in a cl ass of post-AGB 
objects, the so-called w ater fountain sources dlmaiet all 12001: 
Bobo ltz & Mar vel 2005). Additionally, interferometric observa- 
tions of the OH maser region of this class of objects revealed 
aspherical morphologies with either equatorial or bi-conical dis- 
tributions d Amiri et al.l 1201 II) . High resolution observations of 
the SiO masers of the post-AGB object W43A indicated a bi- 
conical outflow dlmai et al.11200 5 1 ). The bi-polar jets observed in 
water fountain sources are related to the onset of asymmetric 
PNe dSahai&Traugerj|1998l) . 

The high fractional linear polarization observed for the SiO 
masers of OH 44.8-2.3 together with the tentative detection of 
circular polarization is potentially an important indication that 
magnetic fields have a significant role in shaping the circum- 
stellar environment of this star. Significant magnetic fields are 
observed in different regions of the circumstellar environment 
of several evolved stars, which indicate the possible role of the 
magnetic field in shaping the CSEs. Observations of H2O masers 
revealed significant field strength for Mira variables and super - 
giants in the range 0.2 G to 4 G (IVlemmings et al.ll2002T, 12005). 
In particular, H2O maser polarimetric observations of the wa- 
ter fountain source W43A showed that the jet is magnetically 
collimated ( Vlemmings et alj|2006l) . Additionally, polarimetric 
observations of OH masers revealed large scale magnetic field 
strength in evolved stars ranging from 0.1 mG to 10 mG (e.g. 
lEtoka & Diamondll2004tlAmiri et al.ll20Toh . 

However, it should be noted that in interpreting the morphol- 
ogy of the maser emission in terms of the underlying physical 
structure a number of aspects have been ignored. This includes 
sufficient column density, velocity coherent path length along 
the amplification path and pumping mechanism. Therefore, it is 
possible that both SiO and OH molecules exist in spherical shells 
but the conditions for the masers mentioned above are not satis- 
fied. 

5.4. SiO emission in OH 44.8-2.3 

The SiO masers of OH 44.8-2.3 exhibit a ring morphology (Fig. 
H]l. Assuming a distance of 1.1 3+0. 34 kpc, the masers are located 
at a distance of ~ 5.4 AU from the surface of the star. The ring 
morphology is similar to the ring patter ns observed previousl y 
for the SiO masers of Mira variables (e.g. Diam ond et alj [l994). 
Moreover, the estimated SiO maser ring radius of OH 44.8-2.3 
is similar to those measured for Mira variables (e.g. 3-7 AU; 
ICotton et ai]|2008l) . However, since OH 44.8-2.3 is a high mass 
loss OH/IR star, the CSE and stellar radius of this star is expected 
to be larger than Mira variables. T he typical stellar radius fo r 
OH/IR stars corresponds to -2.8 AU dHerman & Habingll9 85a). 
This implies that the SiO masers of OH 44.8-2.3 occur at ~1.9 
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stellar radius. This value is similar to the lower end location o f 
SiO masers in Mira variables (2-6 stellar radius: lElitzur| [l992). 
Therefore, even though OH 44.8-2.3 is expected to be larger, the 
SiO masers occur at the same distance from the stellar photo- 
sphere as Mira variables. 

However, the main uncertainty for the discussion above is the 
distance to the star. From the 1612 MHz OH maser observations 
of this source using t he phase lag method a distanc e of 1 .13+0.34 
kpc was measured dVan Langevelde et al.l Il990l) . This method 
assumes spherical symmetry of the OH maser shell. However, 
the distribution of the OH masers of this star (Fig. [2]) does not 
indicate spherical expansion of the OH maser shell. Therefore, 
accurate determination of the distance of this source using the 
parallax method for the SiO maser features of OH 44.8-2.3 are 
necessary to understand the size and location of the SiO maser 
region. 

6. Conclusions 

Our observations indicate a ring morphology for the SiO maser 
region of the OH/I R star OH 44.8-2.3. Assu ming a distance of 
1.13+0.34 kpc (I Van Langevelde et alJ [l990). the masers are lo- 
cated at ~ 5.4 AU from the central star. The ring pattern is simi- 
lar to that observed previously for Mira variables. 

The linear polarization morphology is consistent with the 
dipole magnetic field morphology in the SiO maser region of this 
star. However, we cannot rule out toroidal or solar type field mor- 
phologies. Polarimetric observations of the OH and H2O maser 
regions of the CSE of this star are required to clarify the mag- 
netic field morphology of this star. 

We report a tentative detection of circular polarization at 
~0.7% for the brightest SiO maser feature in the modest spec- 
tral resolution data set. However, due to the increased noise in 
high spectral resolution data set we can not confirm the detec- 
tion. Further polarimetric VLBI observations of the SiO masers 
of this star with more integration time are necessary to clarify 
this. We note that based on polarization studies we can not dis- 
tinguish between the Zeeman and non-Zeeman effects from the 
observations. 

The SiO maser features of OH 44.8-2.3 exhibit two opposite 
arcs. Furthermore, the 1612 MHz OH masers of this star indi- 
cate an elongated shell morphology in the direction where there 
is a gap in the SiO maser emission. Additionally, the direction 
of the magnetic field is parallel or perpendicular to the location 
of the gaps in the SiO maser ring and the OH maser extent. This 
could be taken as a clue that there is a large scale magnetic field 
that imposes a preferred direction on the outflow over scales that 
span 2 orders of magnitude. Because of the timescales involved 
in forming the OH shell, one would then conclude that the mag- 
netic field is important imposing an asymmetric signature on the 
neutral outflow in the OH/IR phase. Furthermore, the high frac- 
tional linear polarization measured for the SiO masers of this star 
could indicate the possible role of the magnetic field in shaping 
the circumstellar environment of this star. However, we can not 
differentiate between the radiative or collisional pumping for the 
SiO maser ring of OH 44.8-2.3. 

Follow up monitoring of the SiO maser region of OH 44.8- 
2.3 can give the distance through parallax measurements and 
this source offers a unique oppo rtunity to check the phase-la g 
distance previously obtained by Van Lang evelde et al.l (jl990). 
Additionally, the future observations will enable us to check the 
evolution of asymmetries in the CSE of this star. 
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